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I Abstract. This article summarizes the structural parameters of the 

Galactic bulge and disk, and discusses the interpretation of the bulge 
q \ microlensing observations and the determination of the Milky Way's lu- 

minous mass from the terminal velocity curve and the Oort limit. The 
bulge is a rotating bar with corotation radius around 4 kpc. The NIR 
disk has a short scale-length. The measured surface density of the local 
! disk is in good agreement with the prediction of a maximum NIR disk 

J> ■ model for the Milky Way. The preliminary new value for the microlens- 

0^ | ing optical depth of clump giant sources is within 1.7a of the prediction 

by the maximum NIR disk model, while the optical depth for all sources 
is still significantly higher. These results imply that cold dark matter 
cannot dominate inside the solar radius. 
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The Milky Way is a normal galaxy seen under a magnifying glass: the ages, 
metallicities, and kinematics of its stars can be observed in much greater detail 
than in any other galaxy. Combined with a model for the large-scale structure 
and potential, dynamical models for specific stellar populations can be con- 
structed which link our local neighbourhood with distant parts of the Galaxy. 
Thus we can study in a unique way the imprints of the evolution and formation 
processes that shaped the Milky Way, and presumably other galaxies as well. 

As is well-known, our viewpoint from within the Galactic disk also has 
some drawbacks, in that many of the Milky Way's large scale properties have 
been difficult to obtain. Only in the last decade has it become possible to map 
the old stellar population in the Galactic bulge and inner disk, using infrared 
observations to penetrate the dust layers in the disk. Models are needed to 
interpret these data and link them to gas kinematical and other observations. 
Two important results are that the inner Galaxy contains a rotating bar /bulge, 
and that the scale-length of the disk is significantly shorter than previously 
assumed. 

Microlensing surveys have now found hundreds of events towards the Galac- 
tic bulge. From such data and from local stellar-kinematic measurements, direct 
constraints on the mass of the Milky Way's bulge and disk can be obtained, 
which are difficult to come by in external galaxies, and which can be used to 
test cosmological models of dark matter halo and galaxy formation. 
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2. Structure of the Galactic bulge and disk 

The best evidence for a rotating bar in the inner Galaxy comes from maps of 
the NIR light distribution, from star counts along various lines-of-sight, from 
the non-circular motions observed in the atomic and molecular gas, and from 
the large optical depth to microlensing towards the bulge. For a more extensive 
review than given here see Gerhard (1999, G99). All length scales given below 
are for a Sun-center distance Rq = 8 kpc. 

Several models for the distribution of old stars in the inner Galaxy have been 
based on the COBE/DIRBE NIR data. These data have complete sky coverage 
but relatively low spatial resolution, they must still be 'cleaned' for residual 
dust absorption, and they contain no distance information, so deprojection is 
not straightforward. The cleaned data show that the bulge is brighter and more 
extended in latitude b at positive longitudes I than at corresponding —I, except 
for a region close to the center where the effect is reversed (Weiland et al. 1994, 
Bissantz et al. 1997). The asymmetry is strongest around |/| ~ 10°. These 
signatures are as expected for a barred bulge with its long axis in the first 
quadrant (Blitz & Spergel 1991). The region of reversed asymmetry at small 
\l\ argues for a bar rather than a lopsided light distribution; see also Sevenster 
(1999). 

Dwek et al. (1995) and Freudenreich (1998) fitted parametric models to the 
cleaned DIRBE data, assuming specific functional forms for the barred bulge and 
excluding low-latitude regions from the fit. In this way the derived outer bulge 
properties are less likely to be influenced by the disk light, but the inner bulge 
and disk profiles are ill-determined. Non-parametric models were constructed 
by Binney, Gerhard & Spergel (1997, BGS), using a Lucy algorithm based on 
the assumption of strict triaxial symmetry, and by Bissantz & Gerhard (2000), 
using a penalized maximum likelihood algorithm. The new model by Bissantz &; 
Gerhard includes a spiral arm model and maximizes eightfold symmetry only for 
the remaining luminosity distribution. In this model the bulge-to-disk ratio in 
NIR luminosity is about 20%, similar to the value given by Kent, Dame k, Fazio 
(1991). Other bar and disk properties from the COBE models are summarized 
below. Physical models for the COBE bar can be found for a range of bar 
orientation angles, 15° ^ (^bar ~ 35°, where ifh&r measures the angle in the 
Galactic plane between the bar's major axis at / > and the Sun-center line. 
<^bar must therefore be determined from other data; see also Zhao (2000). 

The bar is also seen in starcount observations in inner Galaxy fields. Stanek 
et al. (1997) analyzed reddening-corrected apparent magnitude distributions of 
clump giant stars in 12 OGLE fields. The small intrinsic luminosity spread 
(~ 0.2-0.3 mag) makes these stars good distance indicators. The peak of the 
distribution is brighter at / > where the line-of-sight passes through the near 
side of the bar. These fields cover only a small fraction of the sky, but fitting 
a parametric model constrains the bar orientation angle as well as the axis 
ratios and density profile. Nikolaev & Weinberg (1997) reanalyzed the IRAS 
variable population in a similar spirit; here the distance information comes from 
the known range of AGB star luminosities. NIR starcounts have also shown 
longitudinal asymmetries due to the bar (Unavane & Gilmore 1998). Lopez- 
Corredoira et al. (1997, 2000) and Hammersley et al. (1999) have modelled 
the Two Micron Survey Starcounts (mostly bright K and M giants) in several 



3 



strips across the bulge. Structural information on the bulge and disk can be 
derived from these data together with a model for the bright-star luminosity 
function. Ongoing work on deeper surveys (ISOGAL, DENIS, 2MASS) will 
provide important new information on the old stellar population in the inner 
Galaxy; a preview of results is given in van Loon (2000). 

Modelling the HI and CO (l,v) diagrams provides information on the grav- 
itational potential of the bar and disk. Several recent gas flow models have 
followed complementary approaches: Englmaier & Gerhard (1999) modelled the 
gas flow in the potential of the rotating COBE bar of Binney, Gerhard & Spergel 
(1997), Fux (1999) determined a best-fitting N-body-SPH model from a se- 
quence of ab initio simulations, and Weiner & Sellwood (1999) considered gas 
flows in a sequence of analytic model potentials. These simulations produce 
(l,v) diagrams with which many features seen in the observed (I, v) diagrams 
may be qualitatively understood, such as the 3 kpc arm, the non-circular ve- 
locities around the end of the bar, the cusped-orbit shock transition and inner 
X2-disk, the molecular ring and the spiral arm tangent locations. However, no 
model as yet provides a satisfactory quantitative account of the Galactic gas 
disk. A fuller discussion can be found in G99. 

The following subsections contain my best summary of the main bar and 
disk parameters from this and other work. 

Corotation radius: This is the most important bar parameter. The gas- 
dynamical simulations of Englmaier & Gerhard (1999) and Fux (1999) agree in 
their interpretation of the 3 kpc arm as one of the lateral arms close to the bar, 
placing it inside corotation. Sevenster (1999) argues that the 3 kpc arm is part 
of an inner ring, which would also place it slightly inside the corotation radius 
i?CR- The main Galactic spiral arms outside i?cR> ° n the other hand, imply 
an upper limit for Rck, but this is more model-dependent. The gas-dynamical 
models thus give a range of 3 kpc < Rqr 4.5 kpc; a corresponding resonance 
diagram is shown in G99. Dehnen (2000) has interpreted features in the local 
stellar velocity distribution as due to the outer Lindblad resonance with the 
bar, resulting in Rqr = 0.55 ± 0.05i?o> near the upper end of the range from 
the gas-dynamical models. While the match to the Hipparchos data appears 
convincing, it is not clear whether and why this way of determining Rqr should 
work, because the bar is too weak to excite the spiral arms between -Rcr an d 
the solar radius, and these spiral arms should dominate the local quadrupole 
moment (Englmaier & Gerhard 1999). 

Bar orientation From the integrated light alone, physically reasonable models 
can be found for 15° < < 35°. Starcount models give values between 

9?bar = 12 ±6° (Lopez-Corredoira et al. 2000) and 20-30° (Nikoalev & Weinberg 
1997, Stanek et al. 1997). The models of Bissantz & Gerhard (2000) for the 
DIRBE L-band data, when additionally constrained by the clump giant apparent 
magnitude distributions of Stanek et al. (1997), give an optimal i^bar = 15° 
(Figure 1), but i^bar — 15 — 25° is within the uncertainties. The gas-dynamical 
models and the orbit analysis of Binney et al. (1991) are also compatible with 
15° ^ (/?bar 55 35°, depending on whether the emphasis is on the peak in the 
terminal velocity curve, the arm morphology, or the magnitude of the non- 
circular motions near the 3 kpc arm. Finally, microlensing observations favour 
^bar ~ 15° (Zhao & Mao 1996, §3). Thus a good working value is c^bar = 20°. 
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Figure 1. Apparent magnitude distribution of clump giant stars in three 
fields observed by Stanek et al. (1997) with superposed scaled model predic- 
tions from Bissantz & Gerhard (2000) for their model with c/?bar = 15°. The 
model curves have been normalized and shifted along the abscissa such that 
the I = 5.5° and I = —4.9° peaks match best the locations of the observed 
peaks (circles and crosses, respectively) . Note that the relative number, 
width, and asymmetry of the observed distributions in these fields and in 
Baade's window (triangles) are matched well. 



Not consistent with this appear to be the bar model of Hammersley et al. (2000), 
which is based on the identification of a region of strong star formation at / = 27° 
with the nearer end of the bar, and the star count results reported in van Loon 
(2000), which place the near end of a 1.4 kpc size bar at negative longitudes. 
Bar length: Models based on the DIRBE NIR maps find the end of the bar 
around Rgc = 3.2±0.3kpc, when c^bar — 20° (Freudenreich 1998, BGS, Bissantz 
&: Gerhard 2000). This is consistent with with the OH/IR stars (Sevenster 1999), 
IRAS variables (Nikolaev & Weinberg 1997), and the range of -Rcr above for 
a fast bar, while other starcount models use exponential or Gaussian density 
distributions with shorter scale-lengths. 

Bar axis ratios: The parametric DIRBE models give axial ratios of about 
10:3-4:3. This is in agreement with the new non-parametric model of Bissantz 
& Gerhard (2000), whereas BGS had found 10:6:4 without taking into account 
the spiral arms. The starcount models give 10:4:3 (Stanek et al. 1997) and 
10:5.4:3.3 (Lopez-Corredoira et al. 2000). Thus there is good overall agreement 
at around 10:4:3. 

Disk scale-length: In the integrated NIR the radial exponential disk scale Rd 
is significantly shorter than in the optical; numerical values are around 2.5 kpc 
(Freudenreich 1998, BGS) or somewhat shorter (2.1 kpc, Bissantz & Gerhard 
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2000). Hammersley et al. (1999) report satisfactory agreement of their NIR 
counts with a model with Rd = 3.5 kpc, but have not tested other values of Rd- 
Lopez-Corredoira et al. (2000) find that Rd = 3.0 kpc is too short to describe 
the NIR TMSS counts well. Earlier starcount models (Robin et al. 1992, Ortiz 
& Lepine 1993) favour a short disk scale, Rd = 2.5 kpc. It is not clear what 
causes the differences between the various starcount models, and between the 
TMSS starcounts and the integrated NIR light. There may be some interplay 
between the disk scale-length, the bulge profile, and the spiral arm luminosity 
distribution. Further work with spatially complete data will be needed to assess 
this. 

Spiral structure: Spiral arms arc difficult to delineate from our viewpoint 
within the disk. Modelling their (I, t>)-traces directly is complicated by non- 
circular motions (Burton 1973). The most reliable information is available for 
the arm tangent directions of the spiral arms in the gas and young stars (see 
review in Englmaier & Gerhard 1999). The arm segments connecting the tangent 
points are more difficult to determine because of uncertain tracer distances (but 
see Dame et al. 1986). The most widely used spiral model is a four-armed 
pattern following Georgelin & Georgelin (1976). Such a pattern also appears 
most consistent with gas-dynamical models for the Milky Way (Englmaier & 
Gerhard 1999, Fux 1999). Much less is known about spiral arms in the old 
Galactic disk stars (we know from external galaxies that these need not follow 
the blue arms). NIR plane profiles do not clearly show all the spiral arm tangents 
seen in the young component (Drimmel & Spergel 2000), and the contribution of 
young supergiant stars to those that are seen has not been evaluated. Drimmel 
& Spergel (2000) conclude that both a two-armed and a sheared (relative to the 
gas and young stars) four-armed pattern are consistent with the DIRBE data. 
The best models of Fux (1999) contain four stellar spiral arms. 

3. Microlensing towards the Galactic bulge 

Several hundred microlensing events have now been observed towards the Galac- 
tic bulge. These observations give information about the integrated mass density 
towards the survey fields as well as about the lens mass distribution. The most 
robust observable is the total optical depth averaged over the observed fields, 
r. Early measurements gave surprisingly high values T-q ~ 2 — 4 (Udalski 
et al. 1994, Alcock et al. 1997), where r_6 = r/10 -6 . For clump giant sources 
only, which do not suffer from blending problems, Alcock et al. (1997) inferred 
t_6 = 3.9^2 x 10~ 6 from 13 events centered on (l,b) = (2.55°, -3.64°). Us- 
ing a difference image analysis (DIA), Alcock et al. (2000a) recently measured 
t_ 6 = 2.43+^J| for all sources from 99 events centered on (/, b) = (2.68°, -3.35°), 
and from this measurement deduced for the same direction r = (3.23±0.5) x 10~ 6 
for bulge sources only. Finally, in a preliminary analysis of 52 clump giant 
sources in 77 Macho fields, Popowski et al. (2000) found a lower r_6 = 2.0 ± 0.4 
centered on (l,b) = (3.9°, -3.8°). 

It has long been known that axisymmetric models predict T-q ~ 1 — 1.2, 
insufficient to explain the quoted optical depths (Kiraga & Paczynski 1994, 
Evans 1994). Models with a nearly end-on bar as described in §2 enhance r 
because of the longer average line-of-sight from lens to source. The maximum 
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effect occurs for (f> ~ arctan(6/a) when Tb ar /' r axi — (sin2<^>) _1 ~ 2 for 4> = 15° 
(Zhao & Mao 1996). In addition, r increases with the mass and the length of 
the bar/bulge. 

Even so, models based on barred mass distributions derived from Milky 
Way observations (§2) typically give r_6 — 1 — 2 (e.g., Zhao, Spergel & Rich 
1995, Stanek et al. 1997, Bissantz et al. 1997), significantly less than most of 
the measured optical depths. The new bar model of Bissantz & Gerhard (2000) 
gives r_6 = 1.2 for all sources at the position of the DIA measurement and 
r_6 = 1.3 for clump giant sources at the centroid position given by Popowski 
et al. (2000). The mass normalization of the disk and bulge in this model is 
calibrated by assuming constant L-band mass-to-light ratio and by matching 
the predicted gas flow velocities in a hydrodynamic simulation to the Galactic 
terminal velocity curve. As Fig. 1 shows, the apparent magnitude distributions 
for clump giant stars predicted by this model agree closely with those measured 
by Stanek et al. (1997). Thus the model gives a good approximation to the 
distribution of microlensing sources. The quoted optical depths are therefore 
hard to change unless one assumes that the mass distribution of the lenses 
differs substantially from that of the sources. 

Notice that the preliminary new MACHO clump giant optical depth is 
within 1.7cr of the prediction of this bar model. If this clump giant optical depth 
is confirmed, this would be an important step in reconciling galactic structure 
and microlensing observations and would enable us to start using more detailed 
microlensing observables as constraints on Galactic models. On the other hand, 
the recent DIA value is still some 3.2<r away from the model prediction. While 
the NIR model prediction could be slightly increased if the mass-to-light ratio 
were not spatially constant, this is only a ~ 20% effect since limited by the 
terminal velocity curve (Bissantz et al. 1997). Recently, Binney, Bissantz & 
Gerhard (2000) have used general arguments to show that an optical depth for 
bulge sources as large as derived by the MACHO collaboration from the DIA 
value is very difficult to reconcile with measurements of the rotation curve and 
local mass density, even for a barred model and independent of whether mass 
follows light. To illustrate this, the extra optical depth required corresponds 
to an additional mass surface density towards the bulge of some 2000 M / pc 2 , 
comparable to that in the model of Bissantz & Gerhard (2000). Is it possible 
that the DIA measurement is still significantly affected by blending? 

Independent of the resolution of this problem, these results have a further 
important implication. For, determining the mass normalization of the Bissantz 
& Gerhard (2000) model from the terminal velocities implicitly assumes a maxi- 
mal disk. Because the predicted microlensing optical depths are if anything low, 
as much microlensing matter is needed within the solar circle as possible. On 
the other hand, little of the halo dark matter causes microlensing (Alcock et al. 
2000b), so we can not afford significantly more dark mass inside the solar circle 
in the form of CDM particles, say, than corresponding to this maximum disk 
model. Thus the microlensing results argue strongly for a maximum disk. 



7 



4. The maximum disk and local surface density of the Milky Way 

Englmaier & Gerhard (1999) computed a number of gas flow models in the 
gravitational potential of the NIR COBE model determined by BGS, assuming 
a constant joint NIR mass-to-light ratio for the bulge and disk. The actual 
value of Yl can be specified a posteriori. For a maximum disk model it is deter- 
mined by fitting the model terminal curves to the observed HI and CO terminal 
velocities. Some of the models included an additional dark component to prevent 
the outer rotation curve from falling, leading to higher circular velocities at Rq. 
It turns out that the fitted value is insensitive to the circular velocity Vo of 
the Sun implied by the model: within this class of models, the mass of the disk 
and bulge are fixed to within ~ 10%. The luminous component in these models 
accounts for the terminal velocities out to |/| ~ 45°, or R ~ 5.5kpc ~ IRd if 
V = 220 km s -1 , and out to R if V = 180 km s" 1 . 

These maximum disk models predict a surface mass density £© = 45 M Q / pc 2 
near the Sun at Ro = 8 kpc for a local circular velocity of vq = 208 — 220 km s _1 . 
For comparison, the local surface density of 'identified matter' is 48 ± 9 M / pc 2 
(Kuijken & Tremaine 1991, Flynn & Fuchs 1994, Holmberg & Flynn 2000). Of 
this about 23 M Q / pc 2 is in gas and brown and white dwarfs, which contribute 
most to the uncertainty. That the observed and predicted surface density ap- 
proximately agree lends support to the conclusion that the Galaxy indeed has 
a near-maximum disk; the combined observational and model uncertainty is 
about 30% in mass. A NIR disk accounting for only 60% of the rotation velocity 
at 2i?£>, as advocated for spiral galaxies by Rix & Courteau (1999), would have 
only S Q ~ 16 M Q / pc 2 . Turned around, if the NIR disk and bulge are given the 
Tl value implied by the local surface density measurement, then they account 
for the observed terminal velocities in the inner Galaxy. Compared to earlier 
analyses, the main difference is the short disk scale-length (2.5 kpc in the model 
of BGS; see also Sackett 1997) - the Sun is well beyond the maximum in the 
rotation curve from only NIR luminous matter. 

When a cored spherical halo is added to the maximum NIR disk, such as 
to make the Galaxy's rotation curve approximately flat at v c = 220 km s _1 , the 
halo core radius comes out R c ~ 15 kpc. This shows that the Milky Way's dark 
halo is not very strongly concentrated. Integrating the surface density of this 
halo between z = ±1.1 kpc gives 1.1 = 16 M Q / pc 2 . Adding this to the surface 
density of the old NIR disk and the thick disk (~ 9M Q / pc 2 , see discussion in 
G99), the total is Snir + ^th + ^h,i.i = 71 M / pc 2 , whereas the measured total 
Si.i = 71 ± 6M Q / pc 2 (Kuijken & Gilmore 1991). This very good agreement is 
clearly better than one expects, given the uncertainties in both numbers. Note 
that for smaller values of the Galaxy's asymptotic rotation velocity, the required 
amount of halo would be reduced; recall that for v c = 180 km s -1 the terminal 
velocity curve can be fitted within the errors without any added halo (at the 
prize of a falling rotation curve) . A more detailed analysis and comparison with 
cosmologically motivated Milky Way halos is clearly needed. 
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